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ABSTRACT 


A problem in designing semiconductor memories is to provide some measure of 
error control without requiring excessive coding overhead or decoding time. In 
LSI and VLSI technology, memories are often organized on a multiple bit (or 
byte) per chip basis. For example, some 256K-bit DRAW'S are organized in 32Kx8 
bit-bytes. Byte oriented codes such as Reed Solomon (RS) codes can provide 
efficient low overhead error control for such memories. However, the standard 
iterative algorithm for decoding RS codes is too slow for these applications. 

In this paper we present some special decoding techniques for extended 
single-and-double-error-correcting RS codes which are capable of high speed oper 
ation. These techniques are designed to find the error locations and the error 
values directly from the syndrome without having to use the iterative algorithm 
to find the error locator polynomial. Two codes are considered: 1) a d^ n = 4 

single-byte-error-correcting (SBEC) , double-byte-error-detecting (DBED) RS 
code; and 2) a d ^ = 6 double-byte-error-correcting (DBEC) , triple-byte-error- 

detecting (TBED) RS code. 
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I. INTRODUCTION 

Error control has long been used to improve the reliability of computer 
memory systems [1] . The most common approach has been to use a variation of the 
Hamming codes such as the single-error-correcting and double-error-detecting 
(SEC-DED) binary codes first introduced by Hsaio [2] . These codes are particu- 
larly effective for correcting and detecting errors in memories with a 1 bit 
per chip organization. In these memories a single chip failure can affect at 
most one bit in a codeword. ,i-_ 

Large scale integration (LSI) and very large scale integration (VLSI) 
memory systems offer significant advantages in size, speed, and weight over 
earlier memory systems. These memories are normally packaged with a multiple 
bit (or byte) per chip organization. For example, some 256K-bit dynamic random 
access memories (DRAM's) are organized in 32K*8 bit-bytes. In this case a ' 
single chip failure can affect several or all of the bits in a byte, thus ex- 
ceeding the error-correcting and detecting capability of SEC-DED codes. 

Several papers have been written recently trying to extend the SEC-DED 
codes to include byte errors [3-9] . In this paper we investigate the use of 
Reed- Solomon (RS) codes for correcting and detecting byte errors in computer 
memories. RS codes are a class of nonbinary codes with symbols in the Galois 
field of 2 m elements (GF(2 m )). These codes are maximum distance separable (MDS) , 
and thus can provide efficient low overhead error control for byte-organized 
memories, since symbol error correction in GF(2 m ) is equivalent to correcting 
an m-bit byte . 

For computer memory applications, decoding must be fast and efficient. A 
typical RS decoding procedure is to first calculate the error syndromes , then 
use the iterative algorithm [10] to form an error locator polynomial , and 
finally to search for the roots of the error locator polynomial, find the error 
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1 . 
values, and make the actual corrections. The calculation of the error locator 

polynomial is a major step in decoding RS codes, and it remains a bottle-neck 
for high speed decoding, since most errors are single errors and checking, for 
multiple errors is time consuming. High-speed decoding can be achieved by using 
the table-lookup method [1]. However, even for moderate code lengths, the im- 
plementation of table-lookup decoding is impractical, since either a large 
amount of storage or very complex logical circuitry is needed. 

In this paper we investigate some special high-speed decoding techniques 

for extended single-and-double-byte-error-correcting RS codes. These techniques 

are designed to locate and correct the errors directly without having to find 

the error locator polynomial. The occurrences of errors are determined by 

directly testing the weight of the syndrome, denoted by w(£). Let E be the 

number of byte errors. In decoding a single-byte-error-correcting (SBEC) , 

double-byte -error-detecting (DBED) RS code with 3 partiy symbols, if w(s) = 1, 

we show that E = 1. If w(s) * 2, then E _> 2. If w(s) = 3, a simple test is 

required to determine if E = 1 or E > 2. Thus decoding can be carried out in 

parallel, which in effect increases the decoding speed. The detailed procedure 

is presented in Section II. Section III contains our main result, where a 

double-byte-error-correcting (DBEC) , triple-byte-error-detecting (TBED) RS code 

decoding technique is described. Double byte error correction is done by form- 

2 . 

ing a quadratic equation x + x + K = 0, the solution of which gives the two 
error byte locations. The constant K can be determined directly from the syn- 
drome. In this equation, only K contains information about the error loca- 
tions. If a short table is used, with two error locations as two entries cor- 
responding to each value of K, the decoding speed can be made even higher. 
Finally, in Section IV, we summarize our results. 
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= 4 SBEC-DBED CODE 


II. DECODING OF A d . 

mm 

In this section we present a method for decoding an extended Reed-Solomon 

(RS) code over GF(2 m ) with minimum distance d . =4. This code can be used to 

mm 

correct any single byte error and simultaneously detect any double byte error. 
Thus it is called an SBEC-DBED code. Fast encoding and decoding can be acheived 
due to some nice features of the code described below. 


The d . =4 Extended RS Code and Its Properties 

mm r 


It has been shown [11] that there exists an (n+3, n) 
code over GF(2 m ) with parity-check matrix given by 


d . =4 extended RS 

mm 


H = 


£ : 3X3 ! Ml 


( 1 ) 


where I_ _ is the 
— o x 6 

3x3 identity matrix. 
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( 2 ) 


a is a primitive element of GF(2 m ), and n <_ 2 m -l. Because d ^ = 4, the code 

can be used for correcting any single byte error and simultaneously detecting 
all double byte errors [1]. 

From (1) and (2) we see that the H matrix has the following important 
properties. 

1) H is in systematic form. Hence G_ - the generator matrix - is also in 
systematic form: 

G = [ H T i I ] , - (3) 

1 I 

T . 

where is the transpose of H^, and where I_ is the n*n, identity matrix. This 

implies that the parity-check symbols and the syndrome can be generated by the 
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same circuit . 


2) The first nonzero element of every column of H is the unit element 
a^=l . (The advantage of this property will be seen later.) 

3) The number of nonzero elements in each row of H is equal. 

Property 3) simplifies the implementation of the encoder and the decoder 

[ 1 ]. 


Error Correction and Error Detection 

Let v = (Vp, v^ , •••, v n+ 2 ) be a code vector that is written into memory. 
Let r_ = (r Q , r^, • ••, r n+ 2 ^ be the corresponding (possibly noisy) vector that 
is read from memory. Because of possible chip failures, r may be different from 
v. The modulo-2 vector sum 

e = r + v= (e Q , e^ •••, e n+2 ), (4) 

where e. ^ 0 for r. f v. and e. = 0 for r. = v. , is called the error pattern. 

i ill ii* r 

When r is read, the decoder computes the syndrome s, 

_s T = r H T = (v+e)H T = (s Q , s 2 ). (5) 

T 

Since v H = the syndrome s , computed from the vector _r , depends only on the 
error pattern e_, and not on the code vector v. 


a) Single byte error correction 


Let s denote the syndrome corresponding to a single byte error. Then 
~-s 

from (5) we have 


s = e.h. = 
— s 1—1 


( 6 ) 


.th 


where e^ is the error value at location i, and hu is the i— column of H, 

0 < i < n+2 . Note that the first nonzero element of every column of H, is the 
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unit element a^, and e^a^ = e^. Therefore the error value e^ is given directly 
by the first nonzero element of the syndrome. 

The problem of locating the error is reduced to finding a column h. of H 
which satisfies (6) (see Chien [12]). This can be done in the following way. 
Check the elements of the syndrome s to see 


1) if s Q / 0, = s 2 = 0, then i = 0, 

2) if s 1 f 0, s Q = s 2 = 0, then i = 1, 

3) if s 2 j- 0, s Q = = 0, then i = 2. 


Otherwise, since 


s 

— s 


e .h. 


l—i 


1 


" s o‘ 

i-3 



a 


S ! 

a 2 (i-3) 


_ S 2 _ 


we have 



for 3 < i < n+2. 


(7) 


and i gives the error location. Define 



+ s s 
0 2 


( 8 ) 


and note that (7) is equivalent to 


u = 0 for s . j- 0 
3 


b) Double byte error detection 


je{0, 1, 2}. 


(9) 


Let denote the syndrome corresponding to a double byte error. Because 
the code is SBEC and DBED, it follows that [1] 


4 ' • 

for any single and double byte errors. Double byte error detection can be 
done in the following way. If ... ... ; . • . 


( 10 ) 
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( 11 ) 


s i = 0, s i ^ 0, s i f 0, where ± 2 , ± 3 e{0 , 1, 2}, 

or if 

S 1 s 2 

s. t 0, for i = 0, 1, 2, and -± f t ( 12 ) 

1 0 5 1 


then two or more byte errors are detected. Note that (12) implies that 
U * S l 2 + S 0 S 2 * °' 


Summarizing the above discussion, we have the following Decoding Scheme for 
the SBEC-DBED Code (see Fig. 1): 


1) If w(s_) = 0, no errors are detected. If w(s) = 1, E = 1, and 
error correction is done in step 2). If w(^) =2, E _> 2, and errors 
are detected. If w(s_) =3, E _> 1, and decoding proceeds in 

step 3) . 

2) If s. / 0, je{0,l,2), then e. = s. and i = j. 

3 11 ! 

2 

3) Set u = s^ + s q S 2 * If u = 0, E = 1, and calculating 
a 1 3 = s n /s- gives the error location i. Set the error 
value e. = s„. If u / 0, E > 2, and errors are detected. 

l 0 _ 

Figure 2 is a block diagram of the SBEC-DBED decoder. 

III. DECODING OF A d . =6 DBEC-TBED CODE 

mm 

In this section we first present a special high speed decoding technique 
for the DBEC-TBED RS code with d . =6. Then a slightly modified technique is 

applied to decoding the extended RS code with two extra information symbols. 

The d . = 6 RS Code and Its Properties 

mm r ■ 

The generator polynomial for the d^^ = 6 RS code is given by 

2 

g(x) = l (x+a 1 ), (13) 

— i=-2 
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where a is a primitive element of GF(2 m ). The parity-check matrix, H^, of the 
code specified by (13) can be written as 


c* 

1 

-2 

a 

C«‘ 2 ) 

• • • 

-2 n_1 
(a ) 


1 

-1 

a 

(«'b 2 

• • • 

-1 n_1 
(a ) 

»2 = 

1 

1 

1 

• • • 

1 


1 

a 

-(a) 2 

• • • 

, ,n-l 

(a) 


1 

2 

a 

(a 2 ) 

• • • 

2 n_1 
(O 


where n <_ 2 m -l. Because the code has d^^ = 6, it is capable of correcting any 
two or fewer byte errors and simultaneously detecting any combination of three 
byte errors [1] . 

When £ = v + e^ is read, the decoder computes the syndrome £, 

1 = £ H 2 = Cv+e^)H 2 T = e H* = (s_ 2> s_ 1 , s Q , s 1 , s 2 ).'- (15) 


The syndrome corresponding to a single byte error is 


s _ = e . a 

-2 i 


-2i 


s . = e . a 
-1 i 


S 0 = e i 


S ! = e i a 


c 2i 

s 2 = e i a ’ » 


(16.1) 

(16.2) 

(16.3) 

(16.4) 

(16.5) 


where e^ is the error value and i is the error location, 0 _< i <_ n-1, and the 
syndrome corresponding to a double byte error is 1 .. 
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(i7. iy 


-2i -2 j 

e.a + e .a J 

i J 

-i -j 

e.a + e .a 

i J 


e . + e . 

i J 

i j 

e.a + e . a 

i J 


2i 2 j 

e. + e .a 
i 3 


(17.2) 

(17.3) 

(17.4) 

(17.5) 


where 0 <_ i < j <_ n-1. 

Before proceeding, we need to prove some properties of the d ^ ^ = 6 RS 
code which will be used later. 


Property 1 . Let = (s_ 2 , s s qj s ^> S 2^ t ^ ie syndrome corresponding to a 

double byte error. Let N denote the number of zero elements in s^. Then 

N _< 2 , (18) 

and the equality holds in only two cases: 

1) s_ 1 = s 2 = 0; 

2) s^ = s 2 = 0. 

Proof: See Appendix A. 

T 

Property 2, Let s^ = (s s s^, s^, s 2 ) . Then 


s 2 s -2 * S 0 2 * 0 ' 

(19.1) 

V-2 + S -1 S 0 * 0 > 

(19.2) 

Vl * S 2 S -I * °- 

(19.3) 


for all double byte errors. 

Proof : This can be obtained directly from property 1. y-.-. 

Decoding Us ing the Quadratic Equation 

In this subsection we show that the well known quadratic equation over 
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GF(2 m ) can be used to decode the d . = 6 RS code. If a is a primitive element 

mm 

of GF(2 m ); then cT 1 .+ 0 , 0 <_ i < j <_ 2 m -2. From (17.1) and (17.3) we 


have 


det 


e . = 

l 


S -2 “ 


■2j 


det 


-2i -2j 

a a 


-2 j 

S -2 + V 

_• 2 ' 

(a 1 + a ^) 


From (17.2) and (17.3) we have 


det 


s o 1 


e. = 

l 


S -1 a 


-1 


-l -j 

a a 


S -1 + V 

a ^ + a 2 


-3 


Therefore 


S -1 + S 0 a ’ :l S -2 * v' 2 ’ 


-1 -J 
a + a 


(a' 1 +a" J ) 


After multiplying both sides by (a 1 + a ^) ^ 0 and simplifying, the above 

equation becomes 

s ^(a 1 +a' J ) + s 2 ala ■ , + s q = 0. (20) 


In the same way, from (17.3)- (17.5) , we can obtain 


s^c^+a 2 ) + s^a 1 ^ + s 2 = 0. 


Now define 


A 2 

= S 0 + s -i s r 


Y 2 = S 2 S -2 + S 0 ' 


T 3 = S l S -2 + S -lV 


T 4 r Vl + V-l* 


( 21 ) 

( 22 . 1 ) 

( 22 . 2 ) 

(22.3) 

(22.4) 
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Solving (20) and (21) for a 1 +a'* and a 1 ^ , we obtain 


,A i j 
b = a +a 


A i j t 4 
c = a a = — 


for ^ 0* Therefore a 1 and are the roots of 


y + by + c = 0. (2^ 

This is the well-known quadratic equation over GF(2 m ). We will see later that 
it plays an important role in decoding. Therefore we call it the "decoding 
equation". Equation (24) can be rewritten as 


by letting 


where 


x + x + K = 0, 


y = xb. 


K = c/b 2 . 


The formula for the roots of the quadratic equation is (-b ±/b -4c)/2. 
Unfortunately, for finite fields of characteristic two, this formula is not 
applicable because the denominator is zero. However, there are several known 
approaches to solving this equation [10,13,14,15]. The method given in [14] 
is probably the best approach, and we present it in Appendix B. 

Decoding the DBEC-TBED Code 

Suppose that a single byte error with error value e^ at location i occurs, 
From (16. 1)- (16 .5) we see that 

Si t 0, for i = -2, -1, 0, 1, 2, (28.1 


s _± _!l !i i 

S - 2 S -1 " S 0 " S 1 " Q * 
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Note that (28.2) is equivalent to y^ = y^ = y^ = 0. That is, whenever a single 
byte error occurs, s^ ^ 0 for i = -2, -1, 0, 1, 2, and y^ = y^ = Y 4 = 0. From 
(16.3) and (16.4) we have 

s 

a 1 = , (29.1) 

0 


e. 

l 


= s. 


(29.2) 


where i gives the error location and e^ is the error value of a single byte 
error. 

If a double byte error occurs, from property 2 and (22.2)- (22.4) we know 
that y^ 0, y^ ± 0, and y^ £ 0. Therefore b and c in (23.1) and (23.2) exist. 
Hence (24) has two roots, a 1 and . In other words, whenever a double byte 
error occurs, its error locations can be found by solving the decoding equation. 

Since a 1 +a' ) f 0, for 0 <: i < j <_ 2 m -2, when a is a primitive element of 
GF(2 m ), (17.3) and (17.4) imply that 


e. 

l 


det 


det 



V * S 1 

a 1 + a j ’ 


(30.1) 


and 


e . 
J 


+ e. . 

l • 


(30.2), 


where e. and e. are the error values at locations i and j of the double 
1 J 

byte error. 

Let s_ denote the syndrome corresponding to a triple byte error. Then [1] 


’‘if 


( 31 ) 


Based on (31) and properties 1 and 2, we see that if more than two elements 
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x 

of the syndrome £ = (s 2 > s^, s Q , s^, s 2 ) equal zero, but at least one of them 
does not equal zero; or if y 2 , y^, y^ are not all equal to zero, but at least 
one of them does equal zero; or if the decoding equation (24) does not have 
roots in GF(2 m ); then at least three byte errors have occurred. 

Decoding Scheme for the DBEC-TBED RS Code (see Figure 3) : 

XT 

Read £, and calculate the syndrome £ = t_ = (s _, s , s , s.. , s 9 ). 

Let w(y r ) and w (y") denote the Hamming weights of yj = (y^, y^, y ) and 
l" = (y 2 > Y 3 , Y 4 ) , respectively. 

1) If w(£) = 0, no errors are detected. If w(£) = 1 or 2, 

E £ 3 errors are detected. If w(£) = 3 or 4, E > 2, and 
decoding proceeds in step 3). If w(s) = 5, E £ 1, and de- 
coding proceeds in step 2). 

s 

2) Compute y . If w(y ) = 0, E = 1, and calculating a 1 = — 

S 0 

gives the error location i. Set the error value e^ = s Q . 

If w(y’ ) ^ 0, E £ 2 errors are detected, and decoding 
proceeds in step 3). 

3) Compute y" . If w(y") = 3, compute K and T 2 (K) . If 
T 2 (K) =0, E = 2, and we must solve (25) to find the roots 

a 1 and a" 5 . Compute = (s^a^ + s^)/ (a 1 +a- 1 ) and e^ = Sq + e^, 
and correct a double byte error with error values e^ and 
e^. at locations i and j, respectively. If w(y") ^ 3, or 
T 2 (K) = 1, E £ 3 errors are detected. 

Figure 4 is a block diagram of the DBEC-TBED decoder. 

Decoding of the Extended Code 

The parity-check matrix H 2 given in (14) can be extended to form a new 
parity-check matrix given by - 
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(32) 

The code specified by H_ is an (n+2 , n-3) d^ = 6 extended RS code, where 
n _< 2 m -l [16,17,18]. If the errors are confined to location 0 through n-1, all 
the previous results apply. 

Now assume that errors occur at location n or n+1. Then we obtain the 
following results. - 

1) If s_ 2 + s_ x = s 0 = s x = s 2 = 0, (33) 

then a single byte error occurred with error value e^ = s 2 at location n. 

2) If s 2 f s_ 2 = S _ 1 = s Q = s x = 0, (34) 


—3 


»2 


1 0 
0 0 
0 0 
0 0 
0 1 


then a single byte error occurred with error value e n+ ^ = s 2 at location n+1. 


3) If s i f 0, for i = -1, 0, 1, 2, 

and 



(35.1) 


(35.2) 


then a double byte error occurred with error values e. = s_ and e = s _ 
1 l 0 n -2 

SqOi 2l at locations i and n, respectively, where i is obtained from a 1 = 
Note that (35.2) implies that y^ = 0, y^ i 0, and y^ = 0. 

4) If s ± f 0, for i = -2, -1, 0, 1, 


+ 



(36.1) 


and 



(36.2) 


i.e., y^ = 0, y 3 = 0, and y^ £ 0, then a double byte error occurred with error , 

values e. = s„ and e = s_ + s n a 2i at locations i and n+1, respectively,, 
l 0 n+1 2 U 

i ^ 1 

where i is obtained from a = — — . ' - • - 

S 0 
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(37) 


5) If s_ 2 f 0, s 2 f 0, and s_ 1 = s Q = s 1 = 0, 


then a double byte error occurred with error values 
locations n and n+1, respectively. 


e 

n 


s - and e , 
-2 n+1 


at 


Now we combine the discussion in this subsection with that of the previous 
subsection to obtain the following. 


Decoding Scheme for the Extended DBEC-TBED Code (see Figure 5): 

T T 

From the vector r, compute the syndrome £ = _r = (s 2 , s s q» s i» s 2^ 

Again let w(s), w(_y' ), and wQy") denote the Hamming weights of s = 

T 

( s _ 2 > s _i> s 0 > s l > s 2 ^> > I' = 7 3> Y 4 ), and y" = (y 2 , y^, y^) respec- 

tively. 

1) If w(s) =0, no errors are detected. If w(s) = 1, E >_ 1, and 

decoding proceeds in step 2). If w(£) = 2, E _> 2, and decoding 

proceeds in step 3). If w(£) = 3, E >_ 2 , and decoding proceeds 

in step 6). If w(s_) = 4, E _> 2, and decoding proceeds in step 5). 

If w(£) =5, E £ 1, and decoding proceeds in step 4). 

2) If s 2 / 0, E = 1, and a single byte error is corrected with error 

value e^ = s 2 at location n. If s 2 / 0, E = 1, and a single byte 

error is corrected with error value e , = s_ at location n+1. 

n+1 2 

Otherwise, E _> 3 errors are detected. 

3) If s 2 0 and s 2 f 0, E = 2, and a double byte error is corrected 

with error values e = s _ and e . = s,, at locations n and n+1, 

n -2 n+1 2 

respectively. Otherwise, E _> 3 errors are detected. 

i S 1 

4) Compute y_' . If w(y' ) = 0, E = 1, and computing a = — gives the 

0 

error location i. Set the error value e^ = s^. If w(y' ) ^ 0, 
go to step 5) . ; 
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5) Compute y.* . If w(j') = 1 , then: 


(i) If Y 3 t 0, E = 2, and a double byte erroT is corrected 

with error values e. = s_ and e = s _ + s 0 cf 2 ^ a t 

x 0 n -2 0 

locations i and n, respectively, where i is given by 



(ii) If y ji 0, E = 2 and a double byte error is corrected 

with error values e. = s^ and e n = s_ + s„a 2 ^ at 

i 0 n +1 2 0 

i S 1 

locations i and n+ 1 , respectively, where a = — . 

S 0 

(iii) If i 0, E _> 3 errors are detected. 

If w(y ) 1, E _> 2 errors are detected, and decoding proceeds in 

step 6) . 

6 ) Compute _y" . If w(y") = 3, compute K and T 2 (K). If T 2 (K) = 0, 

i i 

E = 2 and we must solve (25) to find the roots a and a and 
then correct a double byte error with error values e^ = 

(s^or 1 + s^)/(a 1 +a J ), e^ = s^ + e^ at locations i and j, respec- 
tively. If w(j") ^ 3, or T 2 (K) = 1, E > 3 errors are detected. 


IV. CONCLUSIONS 

We have presented new decoding techniques for two byte oriented RS codes. 
These decoding techniques are based directly on the syndrome, and do not in- 
volve applying the iterative algorithm to find the error locator polynomial. 
Hence high-speed decoding can be achieved, making these codes well suited for 
correction and detection in byte-organized computer memory ssytems such as LSI 
and VLSI chips . 

The'd^^ = 4 code is capable of single-byte-error-correction (SBEC) and 
double-byte-error-detection (DBED) and can be extended to include three addi- 
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Figure 2. Block diagram of a SBEC-DBED decoder 
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Figure 5. Extended DBEC-TBED decoder error locations calculator 















APPENDIX A 


Proof of Property 1 : It can easily been seen that the vectors (a -2 '*', cT 2 * 1 ), 

O i , a -*), (1, 1), (a 1 , a- 1 ) and (a 2i , a 2 - 5 ), where 0 <_ i < j <_ 2 m -2, are always 
pairwise linearly independent except for the following two pairs: 

1) (a"\ a' j ), (a 21 , a 21 ); 

2) (a 1 , a J ), (a" 21 , a' 2 ^. 

These two pairs are linearly independent for some values of i and j . 

First we show that if s^ = 0, then ^ 0, k = -2, -1, 1, 2. Suppose s^ = 0 
for some k ^ 0. From (17 . 1)- (17 .5) , we have 


1 

a 

if) 


i 

o 


“ 1 


i — 

r-H 

L 

_ s k_ 

ii 

I 

o 

= e. 

l 

ki 

a 

+ e . 
3 

1 

P 

t_l. 

1 


lei lei 

where e^ / 0, e^ / 0 , and k = -2, -1, 1, 2. But (1, 1) and (a , a J ) are lin- 
early independent, and this implies that the above equation is impossible. Hence 
s k t 0, k = -2, -1, 1, 2. 

Next we show that if s ^ = 0 (or s^ = 0), then s^ ^ 0, k = -2, 0, 1, and 
s 2 (or s ^) can be either zero or nonzero. It is easy to show that s^ ^'0, 
k = -2, 0, 1, in the same way as above. Because (a ^ , a ^ ) and (a 21 , a 2 - 1 ) are 
linearly dependent for some i and j, there exists B-^ ^ 0, 3 2 i 0, B-^, B 2 eGF(2 m ), 
and some i < j , such that 


0 


-i 

a 


i 

p 

i 

» — i. 

i 

_ 0 _ 

- e l 

+2i 

a 

+ 6 2 

+2 j 
a 


Let e^ = B^ and e^ = F rom (17.2) and (17.5) we see that the above equation 

becomes . br-c _s.-» 


A-l 



Therefore, s ^ = 0 for some i and j. 

By exactly the same argument as above, we can prove that if s^ (or s ) = 0, 
then s k 5 * 0, k = -1, 0, 2, and s ^ ( or s ^) can be either zero or nonzero. This 
completes the proof that N = 2. Q.E.D. 
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APPENDIX B 


In this appendix we present a method for solving the quadratic equation (25) 
which is based on [14]. 

Let g be any element of GF(2 m ), and define 

A m_1 2 i 

T,(B) * I B . (B.l) 

i=0 


T 2 (S) is known as the trace of $. It is either zero or one [14]. Equation (25) 
has solutions in GF(2 m ) if and only if T 2 (K) = 0 [10,15]. For even m, define 

T. (8) = I B , m even. (B.2) 

i=0 


If (25) has solutions, T^(K) is either zero or one [14]. 

Suppose T 2 (K) = 0, i.e., (25) has solutions. Let x 1 be a solution of (25). 
Then x 2 = 1 + x^ is the other solution, and we have the following results [14]: 

1) m odd 

9 3 9 i 

x x = Ik 2 = Ik 2 

jcJ id 


where I = {1,3,5, •••, m-2}, J = {0,2,4, •••, m-1}. 


2) m E 2 modulo 4 

(m-6)/4 ? ~2+4i 

x = l (K+K ) Z , for T (K) = 0, (B.4.1) 

i=0 

(m-6)/4 7 _2+4i 

X 1 = a l + J ( K+K ) , for T 4 (K) = 1, (B.4.2) 


where is a solution of the equation + 1 = 0. 


3) m e 0 modulo 4 

_ _m-l (m/4)-l ? 2i+m/2 

x = S + S + K (1 + l K ) , for T (K) = 1, (B.5) 

J" - ' Z A ^ 
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where 


(m/ 4) — 1 (m/4)-l r7 2i-l+m/2 2j-2. 

S = l l K (2 + 2 

j=l i=j 


,m. 


For T 4 (K) = 0, select an element g of GF(2 ) such that T (g)- = 1 } compute 
2 2 

K 1 = g+g , and solve z +z+K 1 +K=0 using (B.5) with K replaced 
by K 1 + K. Then x^ = g+z^ is a solution of (25), where z^ is obtained 
from (B.5). For m = 4, 8, 12, (B.5) reduces to the following forms: 


. _.8 „12 
m = 4, x^=K +K ; 


0 V 3Z „66 „129 132 

m=8, x^=K +K + K +K ; 


to „2048,.. „64 „256 ..1024, 

m=12, x^=K(l+K +K +K ) 


+ K 129 + K 258 ♦ K 513 + K 1026 + K 516 + K 1032 , 
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OPTIMAL SHORTENED BCH CODES FOR COMPUTER MEMORY SYSTEMS 


Abstract : This paper presents a method for constructing optimal and nearly 

optimal shortened BCH codes which are suitable for applications to computer 
memory systems. The optimal codes we found minimize the hardware required for 
implementation compared with all the other shortened BCH codes. 

I. INTRODUCTION 

Error-correcting codes are widely used to improve the system-level relia- 
bility of computer main storage or control storage. The IBM system 7030, built 
in 1961, was the first IBM computer system to use a single-error-correcting 
and double-error-detecting (SEC-DED) Hamming code with minimum distance d m ^ n = 4 
for its core memory [1], However, core memories are very reliable, especially 
since the technology has advanced to a mature state. In the 1970's, semicon- 
ductor memories were used to replace core memories. Semiconductor memories are 
faster than core memories in speed; however, they are less reliable than core 
memories due to their high density per chip and their exposure to radiation, 
which induces soft failures (errors). As a result, the use of error-correcting 
codes for improving semiconductor memory reliability became a standard design 
feature. The improvement in reliability is especially evident when the memory 
system is organized on 1-bit-per-card basis (bit-oriented memory) . With this 
organization, most error patterns (or multiple-bit failures caused by a mal- 
function) on each card appear as if they were single errors. 

The most commonly used error-correcting codes have been the minimum odd- 
weight-column SEC-DED codes first constructed by Hsiao [2]. Hsiao's construc- 
tion is optimal in the sense that its parity-check matrix, denoted by H^, 
satifies the following requirements. 



1. The total number of l's in the matrix is a minimum. 

2. The number of l's in each row of is equal (or at least close) 
to the average number (i.e., the total number of l's in divided 
by the number of rows in H^) . 

Let be the number of l's in the ith row of the parity-check matrix. 

Let lP 3 ^ be the number Of logic levels required to generate the ith parity- 

(s') 

check bit with b-input modulo-2 adders (or b-input XOR gates). Let L> . be the 
number of logic levels required to generate the ith syndrome bit with b-input 
modulo- 2 adders. Then we have [1,2] 


^ = [log b (N i -l)1 

CD 

i s) - (1 °<W 

(2) 


where TXl denotes the smallest integer greater than or equal to X. From (1) 
and (2) , we see that the two requirements mentioned above minimize the number 
of logic levels required to generate the parity and syndrome bits, and hence 
hardware required for implementation of the code. 

In practice, some computer memory systems require higher reliability than 
that the SEC-DED codes can provide. These codes should have minimum distance 
dmin >4, and should also satisfy the above two requirements, so that fast en- 
coding and decoding, which are the most critical on-line processes in the 
memory operations, can be achieved. Unfortunately, general methods for con- 
structing such codes are still unknown. In this paper we present an algorithm 
for shortening the BCH codes [1]. In section II, we give some lower bounds on 
the total number of l's in a parity-check matrix H, and on : the average number" 


- 2 - 



of l's in each row of H. In Section III, we present our algorithm, along with 
the shortened BCH codes we found. Finally, we summarize our results in Section 
IV. 


II. LOWER BOUNDS ON THE NUMBER OF l'S IN THE PARITY-CHECK MATRIX 

Throughout this paper, we only consider systematic linear codes, because 
they are the most commonly used codes in practice. The code construction is 
best described in terms of the parity-check matrix H. An (n,k) systematic 
linear code is uniquely specified by its parity-check matrix , which is given 
by 

H = [ 4_ k |Q ], (3) 

where 1^ ^ is the (n-k)x(n-k) identity matrix, 

3. = [ V V v ; 3 ( 4 ) 

is an (n-k)xk matrix, and c^, 1 <i <k, is the ith column of Q. 

TH : Let T and A denote the total number of l's in H and the average number of 

l's in each row of H, respectively. Let m be the unique positive integer such 
that 


l t ": k ) 

i=d ‘ -i i 


+ m = k. 


( 5 ) 


man 


where £ satisfies 


£+1 


i Jv’Jl 1 

■l c n : k ) < k < I c ": k ) 

i=d . -1 1 i=d . -1 1 1 


( 6 ) 


mm 


mm 


and where d . is the minimum distance of the code. Then T and A can be lower 
mm 


bounded by 
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T ^min -1 ) ♦ C »-k ), if k < ^ ) . 

min 


l i C n ; k ) + C * + l )m + ( n-k). 


i=d . -1 
mm 


if ■ f ( n - k ) <k i T c n : k ) . 

i=d . -1 i=d . -1 

mm mm 


( 7 ) 


k( W 1 > n . k 

- ' ' tn -k) — *’■ if k i <■ d . -1 > . 
v mm 


m 


i i ( n : k ) + (A+l )B 

i=d . -1 1 

mm , 

5Ho +1 - 


if 


t . »*1 . 

I c n : k ) <k < l c n : k ) . 

i=d ■ -1 i=d . -1 

mm mm 


( 8 ) 


Proof: Since the code has minimum distance d . , the number of l's in each 

mm 


column of Q must be at least d . -1. If k < f , , the total number of l’s 

— mm — k d . -1 ' 

mm 

in Q is at least k( d ^ -1 )• Adding the (n-k) l’s in the identity part of (3) 
to k( ), we obtain the first part of (7). The second part can be proved 

in a similar way. (8) follows from the fact that H has (n-k) rows. 

An (n,k) systematic linear code is said to be optimal, for a given d ^ , 
if the total number of l’s in its parity-check matrix meet the lower bound of 
(7), and the maximum number of l’s in each row is TA] , where A meets the lower 
bound of (8). As mentioned in section I, no general method is known for con- 
structing optimal codes with d^ n >4. Instead, we seek a method for shorten- 
ing the BCH codes in’ the hope that shortened BCH codes can satisfy -the two re-' ::r 
quirements of section i I. ; ctd 
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Let 


“l =tl n . k ! 1 


( 9 ) 


be the parity-check matrix of an (n,k) systematic BCH code with minimum distance 


d • , where 
mm 


Q : =t i! C1) > a 2 (1) 


Hk (1) 3 


( 10 ) 


is an (n-k)xk matrix. The shortened BCH code is an (n-w,k-w) code specified by 
its parity-check matrix 

— w+1 -m-k : ^w+1 3» ( n ) 

where 

(w+1) 


Q = r (w+1) (w+1) 

*w+l L SLl 1 H-2 


’^k-w 


] 


( 12 ) 


is an (n-k)x(k-w) matrix which is obtained by deleting w columns from Q. , where 

0<w<k. The shortened code has minimum distance at least d . . 

mm 


Define 


= (cw^, cw 2 (1) , . ..,cw k ^3) 

as a column weight vector , where is the Hamming weight of the ith column 


of Q^. Arrange ^ , 1 <_ i <k, in such a way 


that 


CIV. ^3 < CW. (1 3 < ... < CiV. 

h " x 2 - - X k 


If the i, -th, i, -th,..., i, th columns, i.e., the w heaviest columns, 
k-w.+ lr — .-.-k“W+2 — k — 

are deleted from Q. , then the total number of l's in H .is minimized. Let 
—1 — w+ 1 

T. and A. denote the total number of l's in H , and the average number of l's 
1 1 —w+1 

in each row of H w+1 , respectively. Then we have 


k-w 

Till CW. (1) + n-k, (13) 

j=l j 
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and 


k-w 


A 1 1 [ l CWi (1) / (n-k) ] +1 . 
3=1 3 


( 14 ) 


A shortened (n-w, k-w) systematic BCH codes is said to be optimal for a given 

d . if T, meets the lower bound of (13) and the maximum number of l's in 
mm 1 

each row of H „ is FA, 1 , where A, meets the lower bound of (14) . 

— w+ 1 11 


III SHORTENED BCH CODE CONSTRUCTION ALGORITHM 
k 

( ) shortened (n-w, k-w) systematic BCH codes can be obtained from an 

(n,k) systematic BCH code. To obtain the best code, the most straightforward 

k 

approach is to try the ( ) possibilities, and pick up the best one. Unfor- . 
tunately, even for moderate k and w, this approach is impossible due to the huge 
amount of computation. In the following, we present an algorithm for construct 
ing shortened BCH codes. Our results show that some optimal and many nearly 
optimal shortened BCH codes can be found easily by this algorithm. 

Define 


RW (1) = ( Rivf 1 ^, RW, 1 ; 1 ? ..., RW ^ ) 

— 1 2 n-k 

as a row weight vector , where Rwj 1 ^ is the Hamming weight of the ith row of H^. 
The algorithm first sets w, the number of columns to be deleted from Q^, and then 
selects an arbitrary nonnegative real number y > 0. 

Step 1 . Calcualte RW 1 ' 1 ' 1 and 'CW^ . 

Set 


RW^ = max { RW. 1 < i < n-k } 
max . 3 — — 


and calculate 


RW^ 

= — ty, u , 1 < i < n-k, 

> RW (1) “ “ 

max 
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e 


CD 


CD 


(i) 


and 


,CD 

n-k [ 


^ C1) =CW £ (1) -g £ (1)T .A (i:) , 1 < £ < k. 


n)T fll 

where g £ is the transpose of q £ , 


If 

♦ J 15 > ^ , 1 1 j 1 k, 

then the &th column of Qj is deleted. By deleting the ilth column of , we 
obtain an (n-k)x(k-l) matrix 


a 2 ■ (a{ 2 h ^ 


sill 3 


and a new (n-k)x(n-l) parity-check matrix 
-2 = ^hi-k : 

In general, after i < w steps, we have an (n-k)x(n-i) parity-check matrix 


where 


— i+1 ^n-k . Si+i^ 


Q. . , (i+1) (i+1) 

*i+i = tai > S2 


y • • • > 


3k-i - 1, 


which is an (n-k)x(k-i) matrix obtained by deleting i columns from 
Step i+1 . Calculate 


RIV 


(i + D _ [riv ( ^ + D nw(i + D 


, RW, 


1 RW^ + ^^1 

’ *’ n-k J> 



and 


cw ( - 1+1 - ) = ccwJ 1+1 - ) , CW^ 1+1 \ 


CW. 


(i+1) 

k-i 


)• 


Set 


RIV^ + ^ = max {RW.^ + ^, 1 < j < n-k}, 
max . i — J — 

J 


arid calulate 


J 


RW. ( - 1+1 ' ) 

_J 


RW 


(i+1) 

max 


• y. 


1 <_ j < n-k. 


* (i+D 

X 1 

e 



* 


X 

e 


Ci+l) 

n-k 


and 


(i*i) 




1 < Z < k-i. 


If 


(i+1) 


I*- 


(i+1) 


1 <_ j <_ k-i. 


then the & th column of is deleted. Form an (n-k)x(n-i-l) parity-check 

matrix 


— i+2 ’ Un-k'i W 


•i+2 • 


where 


0 = rq (i+2) q (i+2) q (i+2 V r. . 

—i+2 L ^1 * *2 ’ ^k-i-l J, ‘ 

is an (n-k)x(k-i-l) matrix obtained by deleting the ftth column from 
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If i+l.:<w, the above procedure is continued. If i+l=w, H.. the desired 

parity-check matrix of an (n-w, k-w) shortened BCH code. A block diagram of 
the algorithm is shown in Fig. 1. 

In this algorithm, p is a key parameter in finding the best shortened BCH 
codes. The larger the value of y, the smaller the maximum row weight , 

but the total number of l's in the parity-check matrix is larger. On the other 
hand, if y=0, the w columns of with the most l's will be deleted, thus mini- 
mizing the total number of l's in the parity-check matrix of the shortened BCH 
code. But the requirement that the number of l's in each row be equal usually 
can't be satisfied. Which requirement is more important depends on the parti- 
cular situation, and the value of y can be adjusted accordingly. 

Table 1 gives a list of the parameters of the shortened BCH codes found 

it 

by the algorithm, where n = n-w, k = k-w. The (47,32) and (81,64) codes are 
seen to be optimal according to the definition from Section II, and the other 
codes are nearly optimal. Fig. 2 - Fig. 5 show the parity-check matrices of the 
shortened BCH codes with minimum distance d m ^ n =5, 6, 7, and 8, respectively. 

Fig. 6 is a block diagram of the decoder. The encoder is identical to the upper 
part of the decoder with the input parity bits deleted. 

IV CONCLUSION 

In this paper we have presented a method of constructing optimal and nearly 
optimal shortened BCH codes. Some lower bounds on the parameters of the parity- 
check matrix were given. An efficient algorithm for shortening the BCH codes 
was obtained. Some of the codes found by the algorithm were optimal in the 
sense of meeting the lower bounds on the parameters of the parity-check matrix. 
The algorithm was computed on a VAX/780 computer. For moderate values of n and w 
running time was only a few seconds . 
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TABLE 1 PARAMETERS OF A LIST OF SHORTENED BCH CODES 


f 

n 

1 

k 

w 

d . 
min 

Total number 
of l's in H 
(T x ) 

Lower bound 
on T 1 

Lower bound 
on T 

Average number 
of l's per row 
(A 1 ) 

Lower bound 
on A^ 

Lower bound 
on A 

Maximum number 
of l's per row 

28 

16 

1 

35 

5 

88 

86 

76 

7.33 

7.17 

6.33 

8 

44 

i2 

19 

. 5 
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140 
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15.58 

11.67 . 

17 

46 

32 

81 
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12.86 

10.14 

14 

78 

64 

49 
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401 
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270 

28.64 

28.14 

19.29 

30 

80 
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175 
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368 
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272 

23.00 

22.94 

17.00 

24 

29 

16 

34 
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. 93 

7.77 
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79 
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48 
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445 

439 

335 

29.67 
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31 
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409 
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337 
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24.06 

19.82 

25 

34 

16 

29 
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6.89 

6.67 

6.33 
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50 

32 

13 
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13.22 

11.67 

17 

85 

64 

42 
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. 405 

27.38 

27.19 

19.29 

29 

88 

64 

167 
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586 

408 

24.63 

24.42 

17.00 

26 

35 

16 

28 
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133 

131 

131 

7.00 

6.89 

6 189 

8 

51 

32 

12 

8 

261 

257 

243 

13.74 

13.53 

12.79 

16 

86 

64 

41 

8 

656 

656 

470 

29.82 

29.82 

21.36 

32 

89 

64 

166 

8 

673 

663 

473 

26.92 

26.52 

18.92 

28 
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Figure 1. Block diagram of the algorithm for 
shortening the BCH codes. 












Fig.. 2.1 Parity-check matrix of a. (28,16) d . =5. code 
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Fig. 2.2 Parity-check matrix of a (44, 32). code 
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Fig. 2.3 Parity-check matrix of a (46,32) d . =5 code 
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Fig. 4.2 Parity-check matrix of a (50, 32) code 
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Fig. 4.4 Parity-check matrix of a £88, 64) d . =7 code 
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Fig. 5.2 Parity-check matrix of a (51, 32) d . =8 code 
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Fig. 5.4 Parity-check matrix of a (89, 64) d . =8 code 
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Figure 6. Block diagram of the decoder. 
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BURST ERROR CORRECTION IN LASER. MEMORY SYSTEMS 
I. INTRODUCTION 

In computer memory systems, errors often occur in bursts. For example, in 
NASA's laser memory system, the disk can contain a double burst error with 
lengths a _< 8 bits and b _< 11 bits or a single burst error with length c _< 14 
bits. In this paper, we show that by interleaving two d ^ . = 6 Reed-Solomon 
(RS) codes with symbols from GF(2 m ), almost all double burst errors of lengths 
a m+1 and b _< 2m+l can be corrected using erasure decoding, as well as any 
single burst error of length c _< 3m+2. 

Code interleaving distributes the error detection and correction burden 
among the two component codes and thus lowers the overall redundancy require- 
ment. Erasure decoding uses the information inherent in the interleaving scheme 
to achieve further improvement in code performance [1] . More importantly, 
higher speed decoding is possible as opposed to using more powerful codes. 

The generator polynomial for the d^^ = 6 RS code is given by [2,3] 

2 

g(x) = l (x+a 1 ), (1) 

i=-2 

where a is a primitive element of GF(2 m ). The parity-check matrix H of the 
code specified by (1) can be written as 



1 



where n <_ 2 m -l is the code length. If the code is used for error correction and 
detection, we say that the decoder operates in the normal mode . The code is 
then capable of correcting any two or fewer errors and simultaneously detecting 
any combination of three errors [2,3]. A high-speed decoding procedure for the 
decoder operating in the normal mode is given in [4] . If, however, the code is 
used for correcting errors as well as erasures, then all combinations of t errors 
and e erasures are correctable provided that [3] 


2t + e < d . , (31 

min’ v J 

and we say that the decoder is operating in the erasure mode. 

*> t 

Let Vq 1 v^ 1 ■ • • v^ be a codeword generated by the i-th encoder. As shown 

in Figure 1, code symbols from each encoder are sent alternately using a multi- 


i i , 12 12 

plexer resulting m the sequence v^ v^ v^ v^ 


1 2 

v ,v , . This sequence is 
n-1 n-1 n 


then read into the memory where burst noise is introduced. 

At the memory output, each code symbol v^ 1 is replaced by its noise cor- 

i 12 12 12 

rupted version r. . The memory output sequence r^ r Q r^ r^ r n j r n j is 

then demultiplexed into the respective codewords for decoding. The two decoders 
can operate in either the normal mode or the erasure mode. When the memory out- 
put sequence is deinterleaved, assume decoder 1 first operates in the normal 
mode, i.e., it does independent decoding. If the errors in r* = (rj ••• r* 
are correctable, then after error correction the error location information is 
fed into decoder 2 through the mode selector, and decoder 2 will operate in the 


erasure mode. That is, it does erasure decoding by using the error location in- 
formation provided by decoder 1. On the other hand, if the errors in r 1 are 
not correctable but can be detected, decoder 1 stops decoding and decoder 2 
starts operating in the normal mode. If the errors in r = (r^ , .. r n _p 

are correctable, then the error location information from decoder 2 is fed into 
decoder 1. By using this information, decoder 1 will do erasure .decoding and 


2 



the errors in r ^ will be corrected. 

To have a better idea of how the decoder works, let us look at some partic- 

12 12 12 

ular cases. Assume that r = (r Q r Q r^ r^ ^r^ ^ contains two bursts of 

lengths a <_ m+1 and b _< 2m+l. Then the burst of length a < m+1 can affect at 
most two symbols in £ and the burst of length b _< 2m+l can affect at most three 
symbols in r, as shown below. 


,12 12 12 
i = ( r o r o r l r l ••• r i r i 


12 12 
r. ,r. _r. r. 

1-1 3-1 3 3 


1 2 i 

r i r i ) 

n-1 n-1' 


in error in error 


The two code words that result from deinterleaving r are shown below. 
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, 1 1 
(r o r l 


1 
r . 
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in error 


1 

• r . 
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in error 


1 

r 

n- 
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2 

r 








in error in error 

Let decoder 1 operate in the normal mode first. Because r^ contains two 

errors at locations i and j, respectively, they can always be corrected. Then 

2 

the symbols at locations i-1, i, j-1, and j in r are erased. Note that the 

symbols at location i, j-1, and j are in error. By operating decoder 2 in the 

2 

erasure mode, the four erasures in £ can be corrected. If, however, decoder 2 

2 

operates in the normal mode first instead of decoder 1, because r contains 
three errors, they will be detected but not be corrected. Then decoder 2 stops 
decoding, and decoder 1 starts operating in the normal mode. The same decoding 
procedure as described above will now follow. k - l _• 


3 



As another example, again suppose that r contains two bursts, but the first 

2 2 2 
burst only affects r. , and the second burst only affects r. , and r. .as shown 

i 1-1 1 

below. 


,12 12 12 
fr r r t •••r.r 
1 0 0 1 1 i i 


12 12 
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Let decoder 1 operate in the normal mode first. Then no errors will be detected 

in r\ and when decoder 2 starts operating in the normal mode the three errors in 

2 > 

_r will be detected but cannot be corrected. Fortunately,. situations like this 

are rare, as we will show in the next -section. 

The decoding procedure for a decoder operating in the normal mode is pre- 
sented in [4]. Therefore, in the following, we focus our attention only on 
erasure decoding. 


II. ERASURE DECODING PROCEDURE 


Let r 1 = v 1 + e 1 be the input to the i-th decoder. 

(i) 


The decoder first com- 


putes the syndrome £ 


i..T , i i.„T i.,T 
r H = (v +e )H = e H 


rq (i) .CD _(i) (i) 

'--2 , s_ 1 , 5> 0 , s 1 , s 2 J. 


(4) 


Without loss of generality, assume that after the deinterleaving of the 

output sequence, decoder 1 operates in the normal mode. Let t and e denote the 

2 

number of errors and the number of erasures in £ , respectively. Then there 
exists the following possibilities: 

(1) No errors are detected in r^. If this is the case, then decoder 2 
starts operating in the normal mode. If t £ 2, the errors will be 
corrected. If, however, t = 3, the decoder will fail. Suppose 
that each error bit within a burst takes on a value of .1 with prob- 
ability 1/2. Then this failure probability is givenrby 1 j p . r. 


4 



p. = 


1 m+1 


m , ra+l-i 

( l (|) 

i=l z 


[1 - (j) ] > 


... 1 i 1 m+l-i 1 m 

{ l [1 - Oy) Ky) Ky) 


2m+l 


m 

I 

i=l 




m 


(m+1) (2m+l) L . “ ^2 


m m+l-i . i 

{ I (y) [1 - Cy) ] >. 


( 5 ) 


(2) One error at location k in r* is detected and corrected. The symbols 

2 

at locations k-1 and k in £ are erased. Then decoder 2 starts 
operating in the erasure mode. Three cases can occur 

(i) e = 2, t = 0. 


(2) £(k-l) £k 

s £ = e k _ ia + e k a 


a = -2, - 1 , 0 , 1 , 2 . 


( 6 ) 


From (6), the error values at locations k-1 and k are given by 


'k-1 


( 2 ) 


( 2 ) 


1 

k-1 


, (2) k (2), . k,- -1. 

= (s£ + a s^ / a (1+a ), 


(7.1) 


and 


( 2 ) 

e k = S 0 + Vl> 


(7.2) 


respectively. 

(ii) e = 2, t = 1, where the error is at location j f k-1, k. 


sj 2) = e k _ 1 a S ' (k " 1 ' 1 + e k a^ k + e..a £j , £ = -2, -1, 0, 1, 2., (8) 


5 



Define 


T 1 = 


"U (2) 

+ a (l+o ) Sq 

2k- 1 

+ a 

s (2) 
-1 ' 

T 2 = 


k n -1, (2) 
+ a (1+a )s£ 

2k- 1 

+ a 

s ( 2 ) 
s o ‘ 

Then 
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error locator 

and the 

error 
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[31: 





« J * T 2 /Tj , 


C9) 


and 


e . 
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( 2 ) 

0 

( 2 ) 

1 

( 2 ) 


J 

2 j 


T 2 a k Ca _1 + 1) 


, k k-1. , j k-1, , ■ j k. 

(a +a J (ot +a ) (a +a ) \ 


, j k-1., , j k, 
(a +a j (a +a ) 


( 10 . 1 ) 


The other error values are given by 


'k-1 


e. + s< 2 > 1 

J 0 

j (2) , 

e . a + s . 1 

_J 1 

1 1 

k-1 k 

a a 


(e. + + eV + sfb / a k (a _1 + 1), 


( 10 . 2 ) 
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and 


( 2 ) 


+ e 


k-1 


+ e. . 
1 


(iii) e = 2, t = 2, where the errors are at locations j and i (^ k- 
respectively. Because 


2t+e = 4+2 = 6 = d . , 

min 


this error pattern in r cannot be corrected, and the prob- 
ability of this failure occurring is given by 


P 2 = 


1 

m+1 


m 


m 


i *■* 

( l [1 - (t) U • 

i=l 


m 


2m+l 


i i (i - do 

i=l 


Kj) 


m+l-i 


1 2 

(3) Two adjacent errors in _r are detected and corrected. Then r con 

tains at most two byte errors, which can always be corrected by 
operating decoder 2 in the normal mode. 

(4) Two nonadjacent errors in r^ at locations i and j are detected and 

2 

corrected. Then the symbols at locations i-1, i, j-1, and j in r 
are erased, and decoder 2 does erasure decoding. 


( 2 ) 


Jt(i-l) Zi £(j-l) Jt j 

' + e. a + e. ,a + e.a , 

l-l i 1-1 J 


J l = -2, -1, 0, 1, 2. 


Solving (12) for the error values yields 
(i) e i _ 1 = A/B, 


where 


. . (2) i+2j-l (2)., j-1 i., j i, j n -1. 

A = (s^'a + )(cr + a ) (a + a )cr (1 + a ) 

(2 ) r 5j-2,, -1. 2 i+3 j , -3. 3i+2j -1., 

+ Sq J [a J (1 + a ) + a (1 + a ) + a (1 t a )] 

(2) r 4j-l ( - 1 -2. i+3j -1. 3i+j 

+ s^ J [a J (1 + a ) + a (1 + a ) + a (1 + a )], 


(10.3) 
1, k), 

(11) 


( 12 ) 

(13.1) 
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and 


B = 0^(1 + a 2 )(a 21 + a 2 ^ ^ + a 1-1 )^ -1 + a 1 ); 

(ii) e ± = C/D, (13.2) 

where 

r r (2) 3j-l f . i-1. 2j . -1. 

C = (Sq + e i _ 1 )a + (s 1 + e i _ 1 a )a J (1 + a ) 

, 2(i-l). j 

+ Cs 2 + e i _ 1 a J )ct J , 

D = (cr*" 1 + a 1 ) (a^ + ct 1 )^ ; 

(iii) e j _ 1 = [Jcs^ + e i _ 1 + e i ) 

+ sp- 1 + a i (e i _ 1 a 1 + e^] / or’ (1 + a 1 ); (13.3) 

(21 

(iv) e. = J + e i _ 1 + e_. + e^ . (13.4) 

(5) Three errors are detected in r 1 . Then decoder 2 starts operating 

in the normal mode, and the error location information is fed into 

decoder 1. Decoder 1 then does erasure decoding, as described in 

(1) - (4) , the only difference being that if an error at location i 
2 1 

in £ is corrected, the symbols in r_ at locations i and i+1 rather 
than at locations i-1 and i are erased. 


Let P be the probability that two bursts of lengths a = m+1 and b = 2m+l 
occur. Then the probability of decoding failure, from (5) and (11), is given 
by 


P E = (m+l) (2m+l) 


m 

t I [1 

i=l 


1 i i m+l-i 

(i) ](i) 


( (i) 


m 


+ 



(14) 
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If m = 8, then (14) yields 


P c = 3.3PxlO . 
b 

Before proceeding, we need to prove the following property. Let e and t 
denote the number of erasures and the number of errors in r 1 , respectively: 

Property : If e = 2 and 

* Ci) 4 s^CsC 1 ) 2 * .»> ,»)) ♦ ,»>(,(*> ,») + ,(i) .»), 

+ s^ 1 ^ (s^ + s^) = 0, 

then t = 0 or t > 1 . 


We have to 

show t jt 1 . 

It 

is sufficient to show that [3] 



s o (l) 

•S’ 

s Ci) 

-2 



= 

•i (i) 


s (i) 

5 -l 

= 0 . 




s i (l) 

s Ci) 

0 


Because 






A™ 

0 1 0 

2 

+ 

s' 1 ' s^>) ♦ s'j 1 ♦ s' 15 




‘ + 

s 2 1 ' (s -l ) * S Q 1) s 2 1>) = 0 


by our assumption, the property is true. 


Based on the above discussion, our decoding procedure can be summarized as 
follows. Assume that is first decoded by decoder 1 in the normal mode. 

1 .2 

(1) No errors are detected in r . Decoder 2 starts decoding r in ■- ’ 


the normal mode. 


9 



• If two or fewer errors occur in r , they are corrected. 

2 

• If three errors occur in r , they are detected but cannot be 
corrected. 


(2) One error at location k in £ is corrected. 

2 

k-1 and k in r are erased. 


The symbols at locations 


( 2 ) 

• If A J =0, two erasures are corrected with error values e‘ , 

k-1 ' 

and e, given by (7.1) and (7.2), respectively. 

K \ 

• If 0, two erasures and one error are corrected with error 

values e^ ^ and e j given by (10.1), (10.2), and (10.3), 
respectively, where the error location j is given by (9). 

(3) Two adjacent errors in r* are corrected. Decoder 2 starts decoding 
2 

r in the normal mode. 


(4) Two nonadjacent errors in £ at locations i and j are corrected. The 

2 

symbols at locations i-1, i, j-1, and j in r are erased. Four 


erasures are corrected with error values e. ,, e. , e. and e. given 

l-l* l j-l j 

by (13.1), (13.2), (13.3), and (13.4), respectively. 


( 5 ) 


Three errors are detected in £ . 
normal mode. 


Decoder 2 starts decoding r in the 


2 1 

• If no errors are detected in £ , stop. The errors in £ cannot be 
corrected. 

2 

• If one error at location k in r is corrected, the symbols at 

1 ( 2 ) 

locations k and k+1 in £ are erased. Replace £ , k-1, and k 

by _s^ , k, and k+1, respectively, in (9) and (10 . 1) - (10 .3) . De- 
coder 1 corrects two erasures and one error with error values e^., , 

e, , and e. given by (10.1) -(10. 3) , where the error location j 
K^*I J 

is given by (9). ; 


10 



If two errors at locations i and j in r are corrected, the 
symbols at locations i, i+1, j, and j+1 in r* are erased. Replace 
i-1 , i, j-1, and j by i, i+1, j, and j+1, respectively, 

in (13. 1)- (13.4) . Decoder 1 corrects four erasures with error 
values e^, e ^ + ^> e . , and e j + ^ gi ven by (10.1) -(10. 4), respec- 
tively. 


11 
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HIGH SPEED DECODING OF EXTENDED 
TRIPLE-BYTE-ERROR-CORRECTING REED-SOLOMON CODES 


I . INTRODUCTION 


Reed-Solomon (RS) codes are a class of nonbinary codes with symbols or bytes 
from the Galois field of 2 m elements (GF(2 m )). They are maximum distance separa- 
ble, and thus can provide efficient low overhead error control for byte -organized 
memories, since symbol error correction in GF(2 m ) is equivalent to correcting an 
m-bit byte. Chen et al. [1] presented a simplified high speed decoding scheme 
for Reed-Solomon codes capable of correcting up to three byte errors in code 
words made up of k data and n-k parity-check bytes. In this paper, we modify 
Chen's scheme to decode extended triple-byte-error-correcting (TBEC) Reed-Solomon 
codes . 

A typical RS decoding procedure is to first claculate the syndrome, then 
find the error location polynomial and search for its roots, and finally compute 
the error values and make the actual corrections. Finding the error location 
polynomial remains the major bottleneck in high-speed decoding of RS codes. Some 
general solutions to this problem are known, such as those described in [2,3]. 

It is also possible to obtain particular solutions for specific applications. 

The method described in [1] is based on the idea of checking for single errors 
first and correcting them prior to checking for multiple errors. Since most 
errors are single errors and checking for multiple errors is time consuming, this 
method is not only high-speed but also simple to implement. 

In practice, it is also desirable that the coding overhead be as low as 
possible without sacrificing the error correcting capability. Extended RS codes 
provide two more information bytes than ordinary RS codes , and therefore a lower, 
coding overhead, while retaining the same error correcting, capability. ; ;>;In tlyls^n 


1 



report, we modify the decoding scheme given in [1] to extended Triple -Byte- Error- 
Correcting (TBEC) RS codes. The decoder first tests if the second syndrome symbol 
Sj = 0. If / 0, the assumption is made that a single byte error has occurred. 
The assumption is verified quickly by the decoder. If the error is not a single 
byte error the decoder goes on to determine if a double byte error has occurred. 

If not, it then goes on to determine if there is a triple byte error. 

II. PRELIMINARIES 

The TBEC Reed-Solomon code has 6 parity-check bytes or symbols , and minimum 

distance d . =7. The code symbols are elements of GF(2 m ), the finite field of 

mm J 

2 m elements, where m is the number of bits in each symbol. The generator poly- 
nomial g(X) of the code is [3] 

5 . 

g (X) = n (X+a 1 ) , (1) 

i=0 

where a is a primitive element of GF(2 m ). The code length n is assumed to be 
less than or equal to 2 m -l. The parity-check matrix of the code can be written 
as 



“ 1 

1 

1 


1 

1 


1 

a 

(a) 2 


r •> n_2 
(a) 

n-1 

(a) 

H = 

1 

2 

a 

2 2 
(a ) 


9 n-2 

CO 

2 n_1 
(a ) 


1 

3 

a 

(a 3 ) 


3 n " 2 
(a ) 

3 n_1 

(a ) 


1 

4 

a 

(a 4 )' 


4 n ' 2 

(a ) 

4 n_1 
(a 4 ) 


_ 1 

5 

a 

(a 5 ) 2 


5 n ' 2 

(a) 

5 n_1 

(a ) 


Let V = (v., v, , •••, v and R = fr n , r,, •••, r y) be the transmitted 
— '•0 1 n-1' — '•0 1 n-1' 

codeword and the received vector, respectively. The difference : between JR. and l 
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( 3 ) 


V is the error vector E_ = (e Q , e.^, • ••, e n _ 1 )> i-e.. 


R = V * E = (v 0+ a 0 , •••, v n _ 1 * Vl ) 


The syndrome vector is, by definition. 


S = R H T = E H T = (S 0 , S 1 , S 2 , S 3 , S 4 , S 5 ), 


( 4 ) 


and depends only on the error vector, not on the particular transmitted code- 
word. From this definition we have 


S. = £e. (or j . 

j 


Let X. 

tion number. 


- J 


or represent the error location j. X_. is called the error loca- 
Knowing X^ is equivalent to knowing the error location j. Then 


5. = Te.X 1 
i j J J 


0 < i < 5 . 


( 5 ) 


The error- location polynomial a(X) is defined as 

a(X) = n(X-X.) = X e + a X e_1 + • • . + a , (6) 

• 1 X G 

1 

where e is the number of errors in the received codeword. The coefficients of 


ct(X) are related to the syndromes by the following equations [2] : 


S. + O- S. - + 
l+e 1 i+e-1 


+ a ..S. . + aS. = 0 
e-1 l+l e l 


0 = i = 5-e. 


( 7 ) 


The extended triple-byte-error-correcting Reed-Solomon code is an 
(n+2, n-4) code, whose parity-check matrix is given by [4,5,6] 




1 0 
0 0 
H 0 0 

0 0 
0 0 
0 1 


( 8 ) 
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where H is given by (2) . We again have 


S R H £ ( r 0 > •••, r n+1 )Hg - C e 0 > e i> *'•> e n+ i^ H E 


C S Q , s 1 , s 2 , s 3 , s 4 , s 5 ) 


It is easily seen from (8) and (5) that 


= 


n-1 

I e . + e 
j =0 3 n 


(9.1) 


S. 

l 


s, = 


n-1 

I e - X 1 
j=0 

n-1 


1 J 


for 1 < i < 4, 


y e .X. + e . . 

j =0 J 3 n+ l 


(9.2) 


(9.3) 


Note that the error at location n only contributes to S^, and the error at 
location n+l only to S^. 


III. THE MODIFIED DECODING SCHEME 

For the extended code, we denote e as the number of errors at locations 
0 through n-1, e' as the number of errors at locations n and n+l, and E = e+e 1 
as the total number of errors in a received vector, respectively. Obviously, 
if the errors are all confined to locations 0 to n-1, (7) is still valid. 

Then we can write out the equations of (7) explicitly for the following three 
cases : 
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Case 1) E = e = 1: 


+ o^S^ = 

5 3 + °1 S 2 = ° J 

5 4 + °1 S 3 = °* 

5 5 + a l S 4 = °’ 


( 10 ) 


} 


Case 2) E = e = 2: 


S 

S 

S 

S 


2 

3 

4 

5 


* °i s i 


* ° 1 S 2 


+ °1 S 3 


* °1 S 4 


+ 


+ 


+ 


+ 


a 2 S 0 


° 2 S 1 


a 2 S 2 


°2 S 3 


0, 

0, 

0, 

0. 


Case 3) E = e = 3: 


5 3 + °1 S 2 + °2 S 1 + °3 S 0 = 

5 4 + °1 S 3 + °2 S 2 + °3 S 1 = °» 

5 5 + °1 S 4 + °2 S 3 + a 3 S 2 = °* 


(ID 


( 12 ) 


Now let us consider the cases when some of the errors are located at posi- 
tions n or n+1. There exist eight possibilities. It follows from (9.1)-(9.3) 
and (7 ) that : 

Case 4) E = e' = 1 at location n: 

5 0 f 0, S ± = 0, for 1 < i <5. (13) 

Case 5) E = e' = 1 at location n+1: . .. . , 

5 1 = 0, for 0<i<4, S 5 ^0. 3 = . . (14) 
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Case 6) E = e’ = 2 at locations n and n+1 : 


S Q t 0, S. = 0, 1 1 i 1 4, S 5 + 0. 


( 15 ) 


Case 7) E = 2 , e = 1 , e' = 1 , at locations j, 0 < j < n-1, and n: 


S 1 + ° 1 S 0 * °* 


S 2 + ° 1 S 1 = °’ 


S 3 + °1 S 2 = °> 


(16) 


S 4 + °1 S 3 = °> 


S_ + a. S. = 0 . 
5 14 


Case 8) E=2,e=l,e'=l, at locations j, 0 < j < n-1, and n+1: 


S 1 + CT 1 S 0 = °* 


S 2 + a l S l = 0) 


S3 + = 0, 


(17) 


S 4 + G 1 S 3 = °’ 


S c +0,3 f 0. 
5 14 


Case 9) E = 3, e = 1, e' =2, at locations j, 0 < j < n-1, n, and n+1: 


S 1 + Ol s 0 f 0, 


S 2 + CT 1 S 1 = °> 


S 3 + °1 S 2 " °» 


( 18 ) 


S 4 + °1 S 3 = °’ 


s 5 + a 1 s 4 ft 0 . 
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Case 10) 


E = 3, e = 2, e' = 1, at locations j , k, 0 < j < k < n-1, and n: 


S 2 + a l S l + °2 S 0 * °» 


S 3 + °1 S 2 + a 2 S l = °* 


S 4 + °1 S 3 + °2 S 2- = °’ 


( 19 ) 


S c + a.S. + a_S_ = 0, 
5 14 2 3 


Case 11) E = 3, e = 2, e' = 1, at locations j, k, 0 < j < k < n-1, and n+1: 


S 2 + a 1 S ] _ + a 2 S Q = 0, 


S 3 + °1 S 2 + °2 S 1 = °* 


S 4 + a l S 3 + °2 S 2 = °‘ 


( 20 ) 


S 5 + °1 S 4 + °2 S 3 ^ °* 


Our problem is to find a simple way to solve for based on the above 
equations. Let S^, and denote the syndrome vectors for single, double, 
and triple byte error patterns, respectively. Then [3] 


S 

— s 


* ^d * -t ' 


( 21 ) 


The following properties of the syndromes are important to our decoding pro- 
cedure. First, suppose that the error locations are confined to positions 0 
through n- 1 . 

2 

Property 1 ; If t 0 and + S 2 =0, then e = 1, or e _> 3. 


Proof : Because ^ 0, e i 0, and we have to show that e ^ 2. From [2], it is 

sufficient to show that the determinant 


A 


A 



+ S 1 S 3 ‘ °' 
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Since S 2 + = 0 by assumption. Property 1 is true. 

2 

Property 2 : If 0, + S 2 =0, and + S 2 S 3 = 0, then e = 1 or e > 3. 

Proof : From Property 1, we only have to show e ^ 3. From [2], it is sufficient 

to show that the determinant 


S 2 2 * * S 1 S 0 

5 3 S 2 S 1 

5 4 S 3 S 2 


" S 0 (S 3 + S 2 S 4 ) + S 1 (S 2 S 3 + S 1 S 4 ) 


+ s 2 cs 2 + = 0. 


2 2 
Since S^S^ + S 2 =0 and S^S^ + S 2 S 3 = 0 by assumption, = S Q CS^ + S 2 S 4 ). 

2 2 

Also, since S^S^ = S 2 and S^S^ = S 2 S^, then S^ = S 2 /S^ and S^ = S^S^/S 2 (unless 

S„= S_ = S . = 0) . Therefore 
2 3 4 ' 

, s 2 s s 

S 3 ’ sf- • “If " S 2 S 4' and 4 3 “ 

2 2 
Property 3: Let A = S 2 + S^S^ f 0, cr^ = (S 2 S^ + S^S^)/A, and a 2 = (S 2 S 4 + S^ )/A 

(a^ and a 2 when e = 2). If = S^ + CT 1^1 + °2^0 = ® or ^2 ^ ^5 + 0 1^4 + C 2^3 = 
then e = 2, or e > 3. 

Proof: Since A / 0, e ^ 0, and e ^ 1 [2], and it is sufficient to show that 

e f 3. Since the determinant 


A 3 = 


S 2 S 1 S 0 

5 3 S 2 S 1 

5 4 S 3 S 2 


" S 0 (S 3 + S 2 S 4 ) + V S 2 S 3 + S lV + S 2 (S 2 2 + S lV 


SqC^A + SjOj^A + S 2 A = A • = 0, 


it follows that e ^ 3 [2]. Similarly, we can show that if D 2 = 0, e f 3. 
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Now let us consider the cases when the errors can occur at any location. 


Keep in mind that an error at location n can only affect and an error at loca- 
tion n+1 can only affect . 

Property 4 : If f 0 and S^ = 0, 1 <_ i _< 5, then E = e' = 1 with error location 

at position n, or E > 3. 

Proof : Since S Q ^ 0, E f 0 and we need only show that E f 2 and E / 3. From 

Case 4), we see that if E = e* =1 with error location at position n, then 4 0 

and S^ = 0, 1 <_ i <_ 5 . E ^ 2 and E j- 3 follow from (21). Similar arguments can 

be used to obtain the next two properties. 

Property 5 : If S^ i 0 and S^ = 0,0_<i_<4, then E = e' =1 with error loca- 

tion at position n+1, or E > 3. 

Property 6 : If S Q f 0, S 5 f 0, and S^ = 0, 1 _< i ^ 4, then E = e' =2 with two 

error locations at positions n and n+1, respectively, or E > 3. 

2 

Property 7 : Assume S^ f 0, S^S^ + S£ =0, and S.^ + S2S^ = 0, and let = 

S^/S^ (o^ when e = 1), = S^ + o^Sq, and = S 5 + o^S^. Then 

(1) If = 0 and = 0, then E=e=l, or E > 3. 

(2) If ji 0 and T ' = 0, then E = 2, e = 1, and e' = 1, with two error 

locations at positions 3 , 0 _< 3 <_ n-1, and n, respectively, or E > 3. 

(3) If = 0 and f 0, then E = 2, e = 1, and e' = 1 with error loca- 
tions at positions j, 0 3 _< n-1, and n+1, respectively, or E > 3. 

(4) If T / 0 and £ 0, then E = 3, e = 1, and e' = 2 with error loca- 

tions at j, 0 < 3 < n-1, n and n+1, respectively, or E > 3. 

Proof : By Property 2, (1), (2), (3), and (4) follow from Case 1), Case 7), 

Case 8), and Case 9)y respectively. 
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2 2 
Property 8 : Let A = S 2 + f 0, = (S 2 S 3 + S^S^/A, a 2 = (S 2 S 4 + )/A 

(a., and a. when e = 2) , D, = S_ + a, S, + a_S n , and D. = S_ + a. S . + a„S_. Then 
1 2 1 2 11 2 0 2 .5 14 23 

(1) If = 0 and = 0 , then E=e=2orE>3. 

(2) If / 0 and = 0, then E = 3, e = 2, and e' = 1 with error loca- 
tions at positions j, k, 0 _< j < k n-1, and n, respectively, or E > 3 

(3) If D = 0 and j- 0, then E = 3, e = 2, and e' = 1 with error loca- 

tions at positions j, k, 0 <_ j < k <_ n-1, and n+1, respectively, or 
E > 3. 

(4) If • D t 0 and j- 0, then E _> e _> 3. 

Proof : By Property 3, (1), (2), and (3) follow from Case 2), Case 10), and 

Case 11), respectively. (4) follows from Property 3, and (1), (2), and (3). 


Property 9 : If f 0 and = 0 for some i and j such that i j- j and 1 <_ i, 

j _< 4 , then E ^ e ^ 2 . . , 

Proof: First note that e t 0. If e = 1, then S. = e^X. 1 ^ 0, for 1 < i <4. 

i 0 3 “ 

Since some S. = 0, e i 1. Hence e > 2 and E > e >2. 

l — — — 

Property 10 : If = S 2 = 0 and f 0, then E _> 3. 

Proof : Property 10 follows from Properties 1 and 9 

Property 11 : If = S 2 = 0 and ^ 0, then E > 3. 

Proof : From Property 10, E 3, and we need only show that E j- 3. It is 

sufficient to show that 
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Since A3 = 0 by assumption, e j- 3 and E ^ 3. 

Property 12 : If = S2 = S3 = 0 and t 0, then E > 3. 

Proof : Since S^, 1 <_ i _< 4, is not affected by errors at locations n and n+1, if 
S^ = S 2 = S^ = 0 and S^ t 0, it is easy to see that equations (10), (11), and (12) 
are not valid. Therefore e j* 1, e ^ 2, and e ^ 3 and we must have E > e > 3. 


Based on the properties presented above, the procedure shown in Fig. 1 can 
be used to calculate the coefficients of the error-location polynomial (if e / 0) 
and to find the error locations (if e' f 0). 


(1) If S^ / 0, then E > e > 1, and the next step is (2). If S^ = 0 and 

S 2 / 0, then E _> e _> 2, and the next step is (3). If S^ = S^ .= 0 and 

S^ i 0, then E _> e _> 3 and the next step is (7). If S^ = S 2 = S3 = 0 

and S^ f 0, then E > 3 (detected errors). If S^ = 0, 1 _< i _< 4, then 

E _> 0 and e = 0, and the next step is (8). 

(2) S^ i 0. Set a = S 2 /S^. If S3 + aS 2 i 0, then E > e > 2 and the next 

step is (3). If S^ + oS 2 = 0 and S^ + 083 f 0, then E e >_ 3 and the 

next step is (5). If S3 + oS 2 = 0 and S^ + 0S3 = 0, then E _> e = 1, 
and the next step is (6) . 

(3) A = S 2 1 2 3 * + S 1 S 3 f 0. Set ^ = (S^ + S^/A, a 2 = (S 2 S 4 + S 3 2 )/A, 

D 1 = S 2 + °1 S 1 + °2 S 0 * Snd °2 = S 5 + °1 S 4 + °2 S 3 * '•> " ■ n 
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If = D 2 = 0, then E = e = 2. 

If / 0 and = 0, then E = 3, e = 2 and e' =1, with one error 
location at position n. 

If =0 and f 0, then E = 3, e = 2, and e' =1, with one error 
location at position n+1. 

If 0 and ^ 0, then E = e = 3 and the next step is (4). 

(4) Set a 3 = D 2 /D 1 , a 2 = o 2 + a 3 (S^ + S q S 3 )/A, and 

°1 = °1 + °3 (S 0 S 2 + S l 2 - )/A * 

(5) Set A = S 2 S 3 + S^S^ and = Sg(S 3 + S 2 S^) + S-^A . 

If = 0, then E > 3 (detected errors). 

If / 0j then E = e = 3, and we set 

a 3 = [(S 3 + S 2 S 4 )S 3 + A S^J/D^ , = (S 2 S 4 + S^S 3 )/A , and 

°2 = C °3 S 0 + °1 S 2 + S 3 )/S l' 

(6) Set = S-^ + aS Q and T 2 = S 3 + aS 4 » 

If = 0 and T 2 = 0, then E = e = 1. 

I 

If / 0 and = 0, then E = 2, e = 1, and e =1, with one error 

location at position n, 

t 

If = 0 and T 2 ^ 0, then E = 2, e = 1, and e =1, with one error 
location at position n+1. 

I 

If ^ 0 and f 0, then E = 3, e = 1, and e = 2, with two error 

locations at positions n and n+1, respectively. 

(7) S x = S 2 = 0, S 3 / 0. 

If Sq = 0, then E > 3 (detected errors). . y 

If S Q £ 0, then E = e = 3, and we set :: 0, then 
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°1 " W °2 m CS 5 + °1 S 4 )/S 3> and °3 = VV 

(8) Sj - S 2 = S 3 = S 4 = 0. 

If Sq = Sg = 0, then E = 0 (no errors). 

I 

If Sq / 0 and S,. = 0, then E = e =1, with the error location at 
position n. 

! 

If Sq = 0 and j- 0, then E = e =1, with the error location at 
position n+1. 

I 

If Sq f 0 and f 0, then E = e = 2, with the two error locations 
at positions n and n+1, respectively. 

An overall diagram for the decoding of extended TBEC Reed-Solomon codes is 
shown in Fig. 2. The symbols read from the storage media are stored in the data 
buffer while decoding proceeds. In the first step of decoding, syndromes S^, S^, 
S2> S^, S^, and S,. are generated. The syndromes are then stored in the syndrome 
buffer. The OR gates at the output of the syndrome buffer determine if errors 
exist. If errors exist, the decoder starts the error correction procedure. The 
error- location polynomial calculator calculates the error-location polynomial 
using the procedure shown in Fig. 1. The error- location calculator 
then finds the roots of the error- location polynomial. This can be done by using 
the methods described in References 7 and 8. Finally the error values calculator 
determines the value of the errors at each error location. The error locations 
and the error values are used to modify the symbols at the output of the data 
buffer . 

IV. CONCLUSION 

In this report we have presented a decoding scheme for extended triple-byte- 
error- correcting Reed-Solomon codes. The scheme is especially, applicable to com- 
puter systems where -high speed decoding is the basic .requirement... - ,• ; . ;r .h - 


13 



REFERENCES 


[1] C.L. Chen and M.Y. Hsiao, "High Speed Decoding of Reed-Soloraon Codes", U.S. 
Patent 4, 142, 174, February 27, 1979. 

[2] W.W. Peterson and E.J. Weldon, Jr., Error-Correcting Codes , second ed. 

Cambridge, MA: MIT Press, 1972. 

[3] S. Lin and D.J. Costello, Jr., Error Control Coding : Fundamentals and 

Applications , Prentice-Hall, New Jersey, 1983. 

[4] T. Kasami, S. Lin, and W.W. Peterson, "Some Results on Weight Distributions 
of BCH Codes", IEEE Trans. Inf . Theory , IT-12 (2), p. 274, April 1966. 


[5] T. Kasami, S. Lin, and W.W. Peterson, "Some Results on Cyclic Codes Which 
Are Invariant Under the Affine Group", Scientific Report AFCRL-66-62, Air 
Force Cambridge Research Labs., Bedford, MA, 1966. 

[6] J.K. Wolf, "Adding Two Information Symbols to Certain Nonbinary BCH Codes 
and Some Applications", Bell Syst . Tech . J. , 48, pp. 2405-2424, 1969. 

[7] C.L. Chen, "Formulas for the Solutions of Quadratic Equations", IEEE Trans . 
Inform . Theory , IT-28, pp. 792-794, September 1982. 

[8] L.R. Welch, "Fast Algorithms for Finding Roots of Cubic Polynomials over 
Finite Fields", Presented at 1985 IEEE International Symposium on Information 
Theory , June 1985, Brighton, England. 


14 



























Figure 2. 


16 








tional information symbols. The d ^ = 6 code is capable of double-byte-error- 

correction (DBEC) and triple-byte-error-detection (TBED) and can be extended to 

include two additional information symbols. The decoding method applies to the 

extended codes with only slight modification. 

Code efficiency is high since only three parity symbols are used in the 

d . =4 code and only five in the d . =6 code. In addition, the basic code 

mxn mm 

length n can be selected to match the organization of the memory (as long as 

< m » 

n = 2 -1) without changing the decoding method. However, efficiency is maxi- 
mized when n = 2 m -l is chosen. 
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